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Abstract 

<N 

The correlations of the QCD Dirac eigenvalues are studied with use of an ex- 
tended chiral random matrix model. The inclusion of spatial dependence which the 
original model lacks enables us to investigate the effects of diffusion modes. We 
get analytical expressions of level correlation functions with non-universal behavior 
caused by diffusion modes which is characterized by Thouless energy. Pion mode is 
Qs, ', shown to be responsible for these diffusion effects when QCD vacuum is considered 

a disordered medium. 



; 1 Introduction 

It is now widely recognized that a chiral random matrix theory (ChRMT) is a use- 
ful theoretical tool to describe non-perturbative aspects of QCD. Like other complicated 
quantum systems0, statistical properties of QCD Dirac operator show high universality 
described by this model. A possible relation between chiral phase transition and Anderson 
localization transition (a metal-insulator transition induced by a random impurity poten- 
tial) Q, which was originally suggested in the context of an instanton liquid model|| H- 
is currently under intensive study. 

In order to discuss such relations, however, the original model is not appropriate be- 
cause the model lacks space-time coordinate dependence and can only describe a 'metallic' 
region. In this region, energy eigenfunctions extend over the whole system. The corre- 
lations of energy levels are very large and can be correctly described by RMT. But, if 
impurity potential becomes strong, the wavefunction tends to localize. The correlations 
become small and the diffusion effects become important. 

Recently, some authors have discussed the limitation of the original model|| [Tj. Com- 
paring with numerical simulations, they found that the deviations from the results of the 
original model could be seen beyond a certain energy scale J7|, §]. Further there were some 
analytical arguments supporting the above results|| |T0[ . 



1 There are some arguments that Anderson transition can only happen in quenched theories (see e.g. 
Rcf.§) 
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This energy scale is termed Thouless energy after a quantity with a similar role in 
Anderson problem. In Anderson problem, Thouless energy E c is defined as D/L 2 where 
D is a diffusion constant and L is a system size. If energy scale we consider is much smaller 
than E c , the spatial structure of a system can be ignored and we can use a random matrix 
model with a proper symmetry as a good approximation to describe statistical properties 
of a system. But once energy scale becomes comparable to or larger than E c , we can no 
longer ignore the spatial dependence of the system. 

The intuitive explanation of above results is as follows; Since a typical distance of a 
diffusing particle with a diffusion constant D during a time period t is \^Di, a time scale 
in which a particle is diffusing throughout a sample with linear dimension L is r c = L 2 /D. 
The corresponding energy scale, l/r c , is Thouless energy. Hence, when we consider times 
much greater than r c (or energies much smaller than E c ), a particle has enough time to 
wander everywhere in the sample and would finally forget from where it starts. Then 
every point in the sample becomes equivalent and the spatial structure of the original 
system can be ignored. 

Besides the above intuitive reasoning, Thouless energy also appears in expressions of 
various statistical quantities in Anderson problem. In the context of impurity perturbation 
technique, Thouless energy is defined as a first excitation energy of fundamental diffusion 
modes, diffuson. If all excitation energies of the modes with finite wave numbers are much 
greater than the energy scale we consider, we can ignore contributions from these modes. 
There remains the lowest mode with wave number which is constant throughout the 
sample. Keeping only the lowest mode amounts to neglecting the spatial dependence. 

Turning to the statistical properties of QCD Dirac operator, the authors of Ref. [7| 
followed the above intuitive argument and suggested that Thouless energy in QCD is 
the right hand side of Eq.(0). Hence if relevant energy scale exceeds this quantity, we 
should expect deviations from ChRMT due to the effects of diffusion modes. At present, 
however, we do not know concrete expressions of these effects, which is necessary to make 
the arguments quantitative. 

In this paper, we aim to go beyond the qualitative argument and to calculate these 
effects of diffusion analytically. Guided by a model used to describe disordered electrons in 
metalsfTT|1, we extend the original chiral random matrix model; a single random matrix is 
replaced by a set of mutually coupled random matrices each of which is localized in a finite 
space-time region. This model is defined in section 2. Partition function is calculated and 
comparison with the nonlinear sigma model is made. 

To calculate the level density in RMT, one usually use the orthogonal polynomial 



methodflT2f or the supersymmetry method|L3], M. In the present paper, we use the 



supersymmetry method because this method is easy to handle the extended model. This 
is formulated in section 3. 

The averaged level density (section 4) and the two-point level correlation function 
(section 5) are calculated perturbatively. In comparison with the non-chiral RMT, we 
confirm that the Thouless energy in QCD is related to the pion decay constant as originally 
suggested. 

In section 6, we also calculate the two-point correlation functions which correspond 
to meson propagators. We find that a pion propagator derived in the present model 
corresponds to a diffuson propagator in Anderson model. Namely, pion mode plays the 
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role of diffusion mode in disordered QCD vacuum. 



2 Chiral random matrix model 
2.1 Original chiral random matrix model 

Original chiral random matrix theory is defined as follows JTJ . Generally, Euclidean Dirac 
operator matrix in gauge field can be written as 



H 



W \ 
J 

U\ — iuj 2 
L)\ + iu>2 



where W is an N x iV complex matrix and Ui is a real matrix. That is, we consider the 
case of unitary ensemble which has the symmetry for fundamental fermions with three 
colors. N goes to infinity in the thermodynamic limit. Matrix W depends on the gauge 
field and treated as random matrix. Gauge field path integral is replaced by Gaussian 
ensemble average |jj as 

■— TfTrWW^) (2) 




'rmt = / aw i i aet ; t exp 



where rrtf are quark masses and £ denotes chiral condensate. This function is expressed 
by the integral of auxiliary field U which is an Nf x Nf unitary matrix[|l| as 

Z RMT = J dU exp (iV£ Re Tr (m q U)) . (3) 

This function is the same as the partition function of nonlinear sigma model without 
space-time dependence, 

Z NLS = J[dU]exp (- J d 4 xCj (4) 

-> J dU exp (VZReTr (m q U)) (5) 

£2 

C = ^Tr (dptf(x)dJJ(x)) - SRe (Trm q U(x)) + ■■■ (6) 

where is a pion decay constant. In the R.H.S. of the arrow, we ignore the space time 
dependence of U(x). That is, ChRMT is a zero dimensional theory. In ChRMT, we 
focus our attention on quasi- zero mode levels of Dirac operator. As long as we consider 
these quasi-zero modes, we can use the partition function of ChRMT without spatial 
dependence. It is known that the microscopic spectral density whicn is the level density 
of quasi-zero modes does not depend on the assumption of Gaussian ensemble average 



0- 



In nonlinear sigma model, the condition for the zero-momentum sector to dominate 
is I/TOtt 3> L [16[] where m n is the pion mass and L is a system size. This condition 
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means Pion compton wavelength is greater than the system size. When l/m n ~ L, 
contributions from nonzero modes are important and we cannot neglect the kinetic term. 
Using Gell-mann-Oakes-Renner relation, this condition is written as follows 

m q ~ L2 • (7) 

The authors of Ref.|| [7j suggested that Thouless energy in QCD is the right hand side 
quantity of Eq.(^). If relevant energy scale exceeds Thouless energy, we should expect 
that the deviations from ChRMT appear, just like the case of Anderson model. This 



diffusion effects are discussed in Anderson model jn]]. We try to do the same thing in 
QCD. 



2.2 Extended chiral random matrix model 



We must introduce the parameter corresponding to the pion decay constant. To consider 
spatial dependence, the random matrix W is replaced by a set of mutually coupled random 
matrices as 

\ 



W 



\ 







vh 







vh 



where W{ is an I x I random matrix and v is a coupling constant Q. Label i means space- 
time coordinate. Coupling matrix v 1; is chosen diagonal for simplicity. 1; denotes an I 
dimensional unit matrix . Eq.(|]) is the expression for the one dimensional system. But we 
can consider any dimensions. The size of one block matrix, I, represents 'internal' degrees 
of freedom. In Anderson model, this size corresponds to an elastic scattering length. 

Janik et al. suggested that this size corresponds to the constituent quark mass||. 
In instanton liquid model, particularly in finite temperatures, the instantons and anti- 
instantons forms molecules. A random matrix model motivated by this molecule formation 
is discussed in Ref.[irS|. The size of these molecules can be considered I in the present 
model. 

A probability distribution of the random matrices is taken to be Gaussian 

P(H) oc exp |^-/S 2 J2 TrWiWf^ 

where S is a parameter which corresponds to the chiral condensate. When n — 1. 
model reduces to the original one. 

Partition function is calculated using auxiliary field as 



(9) 
this 



J RMT 



A'/' 



H dWi J] det 

i f=l 



rrif 




exp -Z£ 2 J2 TrWiW, 



Y[ dAi exp 



=i 



-/S 2 ^Tr^4 



=i 



det 



2 A similar chiral random hopping model was used 
superconductor/normal- metal systems. 




(10) 



to describe 
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where 



A + rrif 



( (Ax + rrif) ■ h 



t l Nfl 







(A n + m f ) ■ h J 
\ 



V o 



(12) 



l N ji / 



and Ai is an Nf x Nf matrix. This integral is evaluated with use of the saddle-point 
approximation. The saddle-point manifold for v — is — Ui/H where Ui is an Nf x Nf 
unitary matrix. We obtain the partition function as follows 



J RMT 



! J dUi exp 



i=i 



ZlJ^TrmqiUi + Uj) 

8=1 

n-1 

+EVZ Tr (C/^t/, + UjU i+1 ) + ■ 



i=l 



(13) 



The coupling parameter v is assumed to be a small quantity and expanded up to the 
second power. This must correspond to the partition function of the nonlinear sigma 
model, 



: NLS 



J\ dUi exp 



Z£Trm ? (^ + C/t) 



i=l 



+& 1 E,Tt(ul 1 u i + uiu w ) + 



i=l 



(14) 



Here, we set the lattice spacing a = 1. Comparing Eq. ([l3|) with Eq. ([l4]) , we find the 
coupling parameter v is related with the pion decay constant as 



J 7T 



(15) 



In the original ChRMT, the pion decay constant is treated as an infinite quantity. We 
note that v is small but is a large number. This means / is a large number. In the 
original model, the matrix size must be taken infinity in order to get the universal results. 
In the present case, the size of one block matrix / must be taken infinity. On the other 
hand, it is not necessary to take an infinite limit of the number of the block n. 



3 Super symmetry method 



To calculate the level density, we use supersymmetry method[]T3|, We use the nota- 
tions and conventions in Ref. [Q. By this method, calculation of ensemble average can 
be performed easily. However at present this method cannot be applied for unquenched 
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calculation due to the presence of the fermion determinant. We use quenched approxima- 
tion for this reason. In original ChRMT, calculation using this method is performed by 
the authors of Ref . [fL9|, f 



21, E2L 23 



Level density can be expressed as 

p(E) = --Imti G {R) {E), 

where the (retarded) Green function is defined as 
1 



G {R \E) 



E+ -H 

H is a 2nl x 2nl matrix. This is derived from the generating function: 
det(E+-H+ J) 



(16) 



(17) 



Z(E,J) 



det(E+ - H - J) 
c^exp \i^(E + — H — Jk)ip 



k 



1 
-1 



where if) is the 4n/-component graded vector, 

^ t = (<p(+) t a(-) t ,x(+) t ,x(-) t )- 



(18) 
(19) 

(20) 



In the above expression, 0(±) is the nZ-component bosonic vector and x(=t) is t ne 
component fermionic vector, 



(±) 



(tf(±),...,#(±)) 



X T (±) = (xt(±),...,X T n(±)) 

= (Xll(±), • • • , Xu(±), Xnl(±), • • • , Xnl(±)) 



(21) 



(22) 



where the symbol ± expresses the chiral structure of the matrix. 

Taking derivatives with respect to the external field J, we can get the Green function. 
Generating function is normalized to unity when J = and it is easy to calculate the 
ensemble average. 

In Eq. ([I]), chiral basis is used. We perform a unitary rotation as a matter of conve- 
nience and get 



H 



W\ —1002 
ZUJ2 —W\ 



(23) 



This basis is used in Ref. [[19 ] . We follow their way and notations. 
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In this basis, Dirac operator matrix is written as 
H = 



jJi + V 


—iuJi 


) 








iu>2 - 


-oji - V 










( OJll 


o ^ 






' ^21 


o \ 
















V 


Uln j 






{ o 


^2n J 


/ vh 







\ 






vh 














'■■ vli 




• 




V o 







/ 







(24) 



V 



where Uu, u>2i is an I x I matrix and Wi = oju + iuj 2 i. 

Performing the ensemble average over the random matrix, Eq.(^) leads to the follow- 
ing expression 



(Z(E, J)} = J dWZ(E, J) exp |^-/S 2 J2 ^WiW^j 



dip exp 



4ZE 



1 n 

Y,^gA 2 (i) +iip\E + -V - Jk)ip 



i=l 



(25) 



Here, A(i) is a 4 x 4 graded matrix, 

A(i) = -^(A z (i)-Z y A z (i)Z y ) 

where ip^ is a 4 dimensional graded vector 
and H z (Tiy) is the 4x4 Pauli matrix a z (g) I2 (cr 



(26) 
(27) 

(28) 



2 ^ s 



Next, we carry out the Hubbard-Stratonovitch transformation. Eq . (P5|) is expressed 
with use of the auxiliary variable Q as 



(Z(E,J)) 



dQdip exp 



i=l 



+i^(E + — V— Jk)ip] . 



£trgQ(i)i(i 



(29) 



A 4 x 4 graded matrix Q is required to have the same symmetry as A. Since A anticom- 
mutes with E y , we must impose the constraint 



{Q(i),Zy}=0. 

As a consequence, tigQ(i)A(i) = \/2 tigQ{i)A z {i). 



(30) 
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Changing the integration variable to ip — > Si/ 2 ?/' and performing the integral over ip, 
the generating function is expressed as follows 



(Z(E,J)) 



dQ exp 



--5>gQ 2 (*) - Z^trgln £+£ 2 - -Q( 



i=l 



i=l 



i=l 

-1 



+lv 2 J2 trg In E + ^ z - -Q(i) E + ^ z - -Q(i + 1) + 



-i 



+ £ trg In £+£ z --Q(f 
i=i v ^ 



(31) 



where t> is treated as a small parameter and expanded in a power series. The external 
field J is a diagonal matrix. 



4 Microscopic spectral density 



4.1 Saddle-point 

We apply a saddle-point approximation to Eq.(|3l|). Since parameter v is a small quantity, 
we consider the saddle-point for v — 0. The saddle-point equation is 



Q(i) + 



1 



0. 



£+££ 2 - Q(z) 

For ~ O(l ), the saddle-point manifold consists only of the single point 



(32) 



^ f-v-(f); 



(33) 



where minus sign in front of the square root is chosen so that it is consistent with the 
sign of the infinitesimal imaginary part in E + . Level density is expressed as follows 



±, mtr ^<z ( £,y)) 



J=0 



2^ 



^ trg kE z Q(i) exp 



i=i 



^EtrgQ 2 (.) + 

i=l 



(34) 



Using the saddle-point approximation, we can get the familiar semicircle law 



where A = -n jlnVL. 

For E ~ (9(/ _1 ), saddle-point equation becomes 

Q 2 « = -l. 



(35) 



(36) 



This is not a single point and we must perform the integration over the saddle-point 
manifold. 
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Level density is expressed as follows 



(P(E)} 



2tt" 



x exp 



-Im 



Q 2 



dQ trg kE z Q(i 



-ZEE ]T trg Q(i)^ z + ZE V £ trg Q(i)Q(i + !) + ■■ 



i=i 



i=l 



(37) 



To proceed further, we must parametrize the saddle-point manifold. In Ref.[RJ, the 4x4 
graded matrix Q is parametrized as 



Q 

T 



T - l (-iZ z )T 

( (1 + t 2 ) 1 / 2 u 

I -a (l + t 2 ) 1 ' 2 



where 2x2 graded matrix t is given by 

t=( a P *) 
\ p ib I 



(38) 
(39) 

(40) 



We omit the indices % in Q(i) for the sake of simplicity of expressions. Here, a and b 
are real bosonic variables and p and p* are fermionic variables. Ref.||19| uses a diagonal 
parametrization of t, 



t = utgu 1 



/ii 

2/i 



exp 



o e 
e o 



In this parametrization, measure <iQ is calculated as 

AdC i 



dQ 



de x de- 



2n 2 (cosh 9\ cos 6 — i sinh 6\ sin 



dadb 



dpdp* 
2vr 



l--(trg t) 2 + - 



where 



sinh — 
2 

sin -. 

2 



(41) 
(42) 

(43) 
(44) 

(45) 
(46) 



In diagonal parametrization, the measure dQ is expressed compactly. In the exact 
calculation for the original ChRMT, this parametrization is used. But we cannot do exact 
calculation in the present case because we must treat all the nonzero modes. Calculation 
becomes very difficult and we must rely on perturbation. The expression of dQ in the 
form of Eq . (^4|) is suitable for perturbative calculation. 
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4.2 Perturbative calculation 

We take the continuum limit of Eq. (j37|) . The continuum limit for a d dimensional system 
is done by the following replacement 



n — > n 



l^l d } V -> I d d x, Q(i + 1) - Q(i) -> ZVQO) 



(47) 



The level density is expressed as 

(p(E)) = -^-Im J [dQ] (J d d xtrg(ki; z Q(x)) 



x exp 



J d d x (s£trgQ(x)£, - lltigQ(x)V 2 Q( 



(48) 



Here, we used the relation v 2 = / 2 /2E 2 / 2 . 

Eq.(|48]) is expressed by 2 x 2 graded matrices t( 

2S -Re J [dt] (J d d x(l+ tig kt 2 (x)) 



(P(E)) 



7T 

x exp 



J d d x (-4iSEtrgt 2 (x) - 2f 2 trgt{x)V 2 t{x) + 



x 1 



1 

w 

1 



d d x(trgt{x)) 2 + 



8P d 



d d xd d y (trg t (x) ) 2 (trg t ( j/) ) 



+ / rfd ztrgt(a;)trgt 3 (2;) + 



(49) 



We perform the perturbative calculation. Keeping the bilinear terms in the exponential, 
the other terms are expanded in a power series of t and we get 



(p(£)> = V ; Re ( ( 1 



7T 



J- | d d xtrg£;£ 2 (x) - ± J d d x{tigt{x)f 



[nl) 



1 



d d xd d ytrgkt 2 (x)(trgt(y)Y + 



8/ M 



2(n/) d / d 

+^ J d d xtigt(x)tigt 3 (x) + 2f 2 J d d xtig (t(x)Vt(x)) 2 + 



d d xcA/(trgt(:r)) 2 (trgt(?/))' 



((■■•))= /[^(•••)exp 



J d d x (-2/ 2 trgt(a;)V 2 t(x) - 4i££trgt 2 ( 



x 



(50) 
• (51) 



The calculation is just a Gaussian integral and can be performed easily with help of the 
following identity for the arbitrary graded matrix A,B: 

((trg At(x)Bt(x))) = U(E;x,y)trgAtrgB (52) 

( (trg At (x) trg Bt (x))) = U(E;x,y)tvgAB, (53) 

where U(E; x, y) is the fundamental propagator of this model, 

IL{E; x,y)=J2 k)e ik(x - y) (54) 



U{E,k) 



Ms* 2 -if) 



(55) 



7m,- . it 
- (n , = 0, ± l,...), A^^pp, 



f 2 
2E' 
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The perturbative calculation is carried out and we obtain the following expression 

In this calculation, we treat the propagator as a small quantity because the pion decay 
constant is a large quantity. This is correct for the nonzero modes but not correct for the 
zero mode because the term Dk 2 /A vanishes at k = 0. As a result, oscillations due to level 
repulsion are not found. This is because we use the perturbation for all modes. Hence, 
for small energy regions, the above expression is no longer valid due to the divergence of 
k = mode at E — 0. In these regions there are general arguments that the microscopic 
spectral density should take the universal form|T0[|. 



5 Two-point level correlation function 

Two-point level correlation function is calculated in the same way. This function can be 
derived from the ensemble average of the product of the generating functions as 

W(E 1 ,E 2 ) = (tTG (R \E 1 )tTG {R \E 2 ) 



'tr (A,) tr (-zL) J W )Z(E 2 , J^)) . (57) 



4 \dJW J \dJ( 2 ) / 

In ChRMT, we have the following identity 

tiG iA \E) = -tiG {R \-E). (58) 

Here, G^ A \E) is the advanced Green function. With use of this identity, the two-point 
level correlation function is expressed as follows 

(piE^piE,)) = -^[W(E U E 2 )+W(-E 1 ,E 2 ) 

+ W{E 1 ,-E 2 )+W(-E 1 ,-E 2 )\. (59) 

Now let us turn to the calculation of W(E\, E 2 ). The generating function is expressed 
with use of graded matrices and vectors as 

Z(E h J {1) )Z(E 2 , J (2) ) = J d^exp [ift(E -H - Jk)^] (60) 

where 

*-($)■ *-(* i)> -(c ;.)• 

As in the previous section, the ensemble average and Hubbard-Stratonovitch transforma- 
tion is performed and we obtain 

I 



Z(Ei, J^)Z(E 2 , J (2 ^)) = J dQexp 



z i=i 



-zf>g ln(££ 2 - 1q(i)) + ^ 2 ]Ttrg(££ z - ±Q(i))-\EZ z - ±Q(i + l))" 1 
i=i ^ i=i ^ ^ 



n A | 
+ ^trg(EZ z --Q(z))- 1 J' il k + 
i=i ^ 



(62) 
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Here, Q is an 8 x 8 graded matrix. For E ~ 0(1 x ), the saddle point is Q(i) 2 = — 1. 
Parametrization is as follows |19fl 

Q = T^i-iZ^Ta (63) 
To = T u T c h- 



- ch 



An 8x8 graded matrix To is separated by the chiral rotation part T c h and the unitary 
rotation part T u as 

/(1 + t 2 ) 1 / 2 ih 

(1+^) 1/2 it 2 
-it x (1+t 2 ) 1 / 2 

o -it 2 o (i + ^)V2 y 

f (l + t 12 t 2 i) 1/2 o 

(l + t 2 iti 2 ) 1/2 
itn 
-ihi o 



V 



T 



\ 





it 21 

[l + t 12 t2i) 1/2 




^12 




[i + t 21 t 12 y/ 2 J 



(64) 



(65) 



Here, 2x2 graded matrix t is defined as 



t 



12 



ai pi 

Pi ih 

c ir\ 

a id 



t 



21 



a 2 pi \ 
P2 ib 2 ) 

c* a* 
—ir]* id* 



(66) 
(67) 



The chiral rotation is due to the chiral structure of the matrix and the unitary rotation 
is due to the calculation of the two-point level density The former is the same as the 
calculation of the one-point level density in section 4 and the latter is the same as the 
calculation of the two-point level density in the standard (non-chiral) RMT. 
Taking derivative with respect to J, we can obtain W(Ei,E 2 ) as 



W(E l7 E 2 ) = ^f Q2= ^QE^S (*:X z ^Q(z))trg (*£ 

l i3 



1 - A 



4 Jq 
x exp 



Q(i) 



-/£][>g£ + £ z Q« +^ 2 X 2 ][>gQ(i)Q(* + l) 



i=i 



i=i 



where 



A = diag(l,-l,l,-l) 



(6* 



(69) 



In this representation, 8x8 graded matrix is aligned as the following block structure: 
chiral, two-point, boson-fermion, 



M 

M(++) 



M(++) M(+-) \ 
M(-+) M(— ) ) 

M( n )(++) M< 12 )(++) 
M ( 21 )(++) M^(++) 



m££(++) Mff 



(ii), 



(70) 
(71) 
(72) 
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We calculate the integral of Eq.(|68|) perturbatively. The measure dQ is separated into 
the chiral part and unitary (two-point) part: dQ = dp(T c h)dp(T u ) . The chiral part is the 
same as the previous section. Jacobian of the unitary part is equal to unity if we take the 
matrix elements of ti 2 ,t 2 \ as the integral variables. Then we get 

dQ = dt x dt 2 (\ - ^trgt? _ higt 2 2 + ■ ■ dt 12 dt 21 . (73) 

The calculation is carried out in the same way as the previous section and yields 
W(E U E 2 ) = -A(nl) 2d i: 2 + 8i: 2 n d (nl) 2d J2(^ 2 (E 1 ,k) + U 2 (E 2 ,k)) 

k 

-i 6 eV) m eW^^, *) + ••■. ( 74 ) 



Two-point level density is expressed as 
(p(E 1 )p(E 2 )) - (p(E 1 )) (p(E 2 )) 

= i £ [ n2 ( E > fc ) + n2 ( -jB > fc ) + n2 ^/ 2 > fc ) + n2 ( _cu /2 ' k ">] + ' ' * > ( 75 ) 

where = (Ei + E 2 )/2 and uu = Ei — E 2 . independent terms are characteristic in 
ChRMT. Comparing this expression with the usual (non-chiral) RMT, we confirm that 
Thouless energy in QCD is 

E c = (76) 



We can calculate number variance from the two-point level density. This is defined as 

rz/2 
-z/2 



T?{z) = [ Z/2 dxdyA 2 {{p{E x )p{E 2 )) - (p(E 1 )) (p(E 2 ))) , (77) 

J -z/2 



where x = E±/A and y = E 2 /A. This is the measure of fluctuations of the level number 
in the finite level interval z. 

Integration is performed and we obtain 

tfW-sEfcfl + ifgjisV (78) 



7T- 



k \ v A 



We note that the ^-dependent terms and a;-dependent terms give the same contributions 
to the number variance. Eq. (|75f) does not treat the zero mode exactly. We expect, 
however, that this approximation does not change the behavior of T?(z) drastically. This 
is because it is known from the results of the usual RMT that what is missing in the 
perturbative result of the two-point level density is the oscillating behavior of the exact 
result and the integration in Eq . (|77|) cancels these oscillations. Taking the continuum 
limit, we obtain 

= h(wf*J M *( 1 + l?) < TO > 
oc z d ' 2 . 
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We find that the number variance depends on the dimensionality of the system. This 
is very different from the original ChRMT where £ 2 (z) ~ logz for large z. The similar 
behavior is find in Ref.[33] for usual (non-chiral) RMT and interpreted as change from 
level correlation of RMT to those of non-correlated energy levels (Poisson statistics : 
Y?(z) ~ z). This is consistent with intuitive discussion. 



6 Two-point correlation function 

The spatial dependence, which we introduced into the original ChRMT in this paper, 
enables us to calculate the correlation between two spatial point. By formal replacement 
of the energy E + to im, the Green function G^ R \E) becomes a quark propagator. The 
ensemble average of a matrix element 



im\ x,y,jj,,v) = \ x,jx 



1 



y,v) (80) 



im — H 

can be calculated from the generating function as 

1 d 

{G(im; x, y, fj,, v)) = -— (Z(im, J)) . (81) 

The external field J is not diagonal in this case. Label x, y represents the position of 
block matrices. Label \i,v represents the 'internal' degrees of freedom which are described 
and by random matrices. We get (G(im; x, y, v)) — if \x — y\ > I. This behavior 
means the size I is a reciprocal of the constituent quark mass. This is consistent with the 
discussion of the Ref.g. 

We can also consider 'meson' propagator in this model. This is the product of the 
Green functions and can be calculated from the generating function as 

(try] G(im; x, y, /i, u)G(im; y, x, v,[i)) = 

1 y ( d d + d d 

d d d d \ 

+ ar , K^Al , )(Z(E,J)) (82) 



and 



dJxy^iui^ "I") dJyxfiu ("I - ) dJxyfiui ) dJyx^yi^ ) 



I tr J2 G(im; x, y, //, is)j 5 G(im; y, x, is, /j)j 5 \ = 
ly^f d d + d d 

d d d d \ 

+ ^ , ^ , + , T i ^ rr, (m J)) . (83) 



dJxyfiu{ ~t~) dJyx^,u{ "I") dJ X yfJ,u(y ) dJy X ^u(y~^~) / 

Here, trace is taken over the chiral degrees of freedom. Summation over fi,u is taken 
because this degrees of freedom is described by random matrices. Eq. (|82"D corresponds to 
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scalar (sigma) propagator and Eq.(p3[) corresponds to pseudo scalar (pion) propagator []. 
We consider the long range correlation between the point i / !/. In the present basis, 
matrix 75 is expressed as 

Eq.(|52"D and ( |83| ) are calculated perturbatively and we obtain 

( tr ^2 G(im; x, y, /i, u)G(im; y, x, /i) ) 
\ fj.,f I 

= 64S 2 / d J d d zU(im;x,y)U( im; y , z)H(im; z , x) + • • • (85) 
( tr G(im; x, y, //, u)j 5 G(im; y, x, u, /i)^ 5 ) = -16l 2d T, 2 Il(im; x, y) H (86) 



where 



U(im;x, y) = U(im, k) exp (ik(x — y)) (87) 

k 

1 



II(irn, k) 



47T (f P + f 
1 

47rg(fc 2 + m2) 



(88) 



Here Il(im; fc) is obtained from II(£'; fc) by the replacement of E by im. The above 
equation shows that H(im; x, y) corresponds to the pion propagator. We find that the 
scalar propagator can be neglected by comparing with the pseudo scalar propagator. 
Pseudo-scalar mode is important. This is consistent with the NG theorem. As was shown 
in Section 2.2, the present model is equivalent to the nonlinear sigma model. Hence it is 
reasonable that the pion propagator appears as a fundamental ingredient in the calculation 
of various correlation functions. The calculation of the two-point correlation function was 
also done in Ref. with use of partially quenched chiral perturbation theory. 

H(E; x, y) is a propagator of the fundamental mode in the perturbative calculation. 
In Anderson problem, this mode is called diffuson. Eq. (|86|) shows the pion mode in QCD 
plays the same role as the diffuson mode in Anderson model. The fundamental mode 
appealing in the chiral perturbation theory, that is, pion, corresponds to the fundamental 
mode appealing in the Anderson model, diffuson. 



7 Conclusions 

In this paper, we examined the level correlations of QCD Dirac operator. We introduced 
the extended chiral random matrix model including the spatial dependence. 

3 We must insert the isospin matrix but this is irrelevant in this calculation. 
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First, we showed the partition function of this random matrix model is equivalent to 
that of a nonlinear sigma model. This equivalence has been previously shown only for 
the O-dimensional case without spatial dependence. As far as we know, this is a first 
explicit derivation of nonlinear sigma model (with derivative terms) from the kind of 
chiral random matrix models. 

To calculate level statistics, we used supersymmetry method developed by Efetov. 
Using the perturbative expansion, we got the expression of the level density and two- 
point level correlation function. 

Comparing the expression of this function with the usual RMT, we confirm that the 
Thouless energy in QCD is E c = f^/2'EL 2 . Number variance was calculated from the two- 
point level correlation function and found to depend on the dimensionality of the system. 
This is very different from the standard universal result and interpreted as change from 
level correlation of RMT to those of non-correlated energy levels. 

For the first time, the correlation functions of pseudo-scalar type as well as scalar 
type between two spatial points could be calculated for the chiral random matrix model. 
The previous work with use of partially quenched chiral perturbation theory^ only gives 
the correlation function of pseudo-scalar type. The reason why scalr type correlation 
function can be calculated in this work is that quark propagators can be defined for the 
chiral random matirx model. 

These functions demonstrate that the pion propagator is nothing but the diffuson 
propagator which is a fundamental mode responsible for diffusion effects in Anderson 
model. Namely, pion mode plays a role of diffuson when we consider QCD vacuum as a 
disordered medium. On the other hand, sigma-like mode does not appear as a fundamental 
mode in this model. 
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